Bacterial Genes in the Aphid Genome: Absence of Functional Gene Transfer from Buchnera to Its Host

Genome reduction is typical of obligate symbionts. In cellular organelles, this reduction partly reflects transfer of ancestral bacterial genes to the host genome, but little is known about gene transfer in other obligate symbioses. Aphids harbor anciently acquired obligate mutualists, Buchnera aphidicola (Gammaproteobacteria), which have highly reduced genomes (420–650 kb), raising the possibility of gene transfer from ancestral Buchnera to the aphid genome. In addition, aphids often harbor other bacteria that also are potential sources of transferred genes. Previous limited sampling of genes expressed in bacteriocytes, the specialized cells that harbor Buchnera, revealed that aphids acquired at least two genes from bacteria. The newly sequenced genome of the pea aphid, Acyrthosiphon pisum, presents the first opportunity for a complete inventory of genes transferred from bacteria to the host genome in the context of an ancient obligate symbiosis. Computational screening of the entire A. pisum genome, followed by phylogenetic and experimental analyses, provided strong support for the transfer of 12 genes or gene fragments from bacteria to the aphid genome: three LD–carboxypeptidases (LdcA1, LdcA2,ψLdcA), five rare lipoprotein As (RlpA1-5), N-acetylmuramoyl-L-alanine amidase (AmiD), 1,4-beta-N-acetylmuramidase (bLys), DNA polymerase III alpha chain (ψDnaE), and ATP synthase delta chain (ψAtpH). Buchnera was the apparent source of two highly truncated pseudogenes (ψDnaE and ψAtpH). Most other transferred genes were closely related to genes from relatives of Wolbachia (Alphaproteobacteria). At least eight of the transferred genes (LdcA1, AmiD, RlpA1-5, bLys) appear to be functional, and expression of seven (LdcA1, AmiD, RlpA1-5) are highly upregulated in bacteriocytes. The LdcAs and RlpAs appear to have been duplicated after transfer. Our results excluded the hypothesis that genome reduction in Buchnera has been accompanied by gene transfer to the host nuclear genome, but suggest that aphids utilize a set of duplicated genes acquired from other bacteria in the context of the Buchnera–aphid mutualism

Aphids acquired symbiotic genes via lateral gene transfer

<p>Abstract</p> <p>Background</p> <p>Aphids possess bacteriocytes, which are cells specifically differentiated to harbour the obligate mutualist <it>Buchnera aphidicola </it>(γ-Proteobacteria). <it>Buchnera </it>has lost many of the genes that appear to be essential for bacterial life. From the bacteriocyte of the pea aphid <it>Acyrthosiphon pisum</it>, we previously identified two clusters of expressed sequence tags that display similarity only to bacterial genes. Southern blot analysis demonstrated that they are encoded in the aphid genome. In this study, in order to assess the possibility of lateral gene transfer, we determined the full-length sequences of these transcripts, and performed detailed structural and phylogenetic analyses. We further examined their expression levels in the bacteriocyte using real-time quantitative RT-PCR.</p> <p>Results</p> <p>Sequence similarity searches demonstrated that these fully sequenced transcripts are significantly similar to the bacterial genes <it>ldcA </it>(product, LD-carboxypeptidase) and <it>rlpA </it>(product, rare lipoprotein A), respectively. <it>Buchnera </it>lacks these genes, whereas many other bacteria, including <it>Escherichia coli</it>, a close relative of <it>Buchnera</it>, possess both <it>ldcA </it>and <it>rlpA</it>. Molecular phylogenetic analysis clearly demonstrated that the aphid <it>ldcA </it>was derived from a rickettsial bacterium closely related to the extant <it>Wolbachia </it>spp. (α-Proteobacteria, Rickettsiales), which are intracellular symbionts of various lineages of arthropods. The evolutionary origin of <it>rlpA </it>was not fully resolved, but it was clearly demonstrated that its double-ψ β-barrel domain is of bacterial origin. Real-time quantitative RT-PCR demonstrated that <it>ldcA </it>and <it>rlpA </it>are expressed 11.6 and 154-fold higher in the bacteriocyte than in the whole body, respectively. LdcA is an enzyme required for recycling murein (peptidoglycan), which is a component of the bacterial cell wall. As <it>Buchnera </it>possesses a cell wall composed of murein but lacks <it>ldcA</it>, a high level of expression of the aphid <it>ldcA </it>in the bacteriocyte may be essential to maintain <it>Buchnera</it>. Although the function of RlpA is not well known, conspicuous up-regulation of the aphid <it>rlpA </it>in the bacteriocyte implies that this gene is also essential for <it>Buchnera</it>.</p> <p>Conclusion</p> <p>In this study, we obtained several lines of evidence indicating that aphids acquired genes from bacteria via lateral gene transfer and that these genes are used to maintain the obligately mutualistic bacterium, <it>Buchnera</it>.</p

Functional Convergence in Reduced Genomes of Bacterial Symbionts Spanning 200 My of Evolution

The main genomic changes in the evolution of host-restricted microbial symbionts are ongoing inactivation and loss of genes combined with rapid sequence evolution and extreme structural stability; these changes reflect high levels of genetic drift due to small population sizes and strict clonality. This genomic erosion includes irreversible loss of genes in many functional categories and can include genes that underlie the nutritional contributions to hosts that are the basis of the symbiotic association. Candidatus Sulcia muelleri is an ancient symbiont of sap-feeding insects and is typically coresident with another bacterial symbiont that varies among host subclades. Previously sequenced Sulcia genomes retain pathways for the same eight essential amino acids, whereas coresident symbionts synthesize the remaining two. Here, we describe a dual symbiotic system consisting of Sulcia and a novel species of Betaproteobacteria, Candidatus Zinderia insecticola, both living in the spittlebug Clastoptera arizonana. This Sulcia has completely lost the pathway for the biosynthesis of tryptophan and, therefore, retains the ability to make only 7 of the 10 essential amino acids. Zinderia has a tiny genome (208 kb) and the most extreme nucleotide base composition (13.5% G + C) reported to date, yet retains the ability to make the remaining three essential amino acids, perfectly complementing capabilities of the coresident Sulcia. Combined with the results from related symbiotic systems with complete genomes, these data demonstrate the critical role that bacterial symbionts play in the host insect’s biology and reveal one outcome following the loss of a critical metabolic activity through genome reduction

Proteinortho: Detection of (Co-)orthologs in large-scale analysis

<p>Abstract</p> <p>Background</p> <p>Orthology analysis is an important part of data analysis in many areas of bioinformatics such as comparative genomics and molecular phylogenetics. The ever-increasing flood of sequence data, and hence the rapidly increasing number of genomes that can be compared simultaneously, calls for efficient software tools as brute-force approaches with quadratic memory requirements become infeasible in practise. The rapid pace at which new data become available, furthermore, makes it desirable to compute genome-wide orthology relations for a given dataset rather than relying on relations listed in databases.</p> <p>Results</p> <p>The program <monospace>Proteinortho</monospace> described here is a stand-alone tool that is geared towards large datasets and makes use of distributed computing techniques when run on multi-core hardware. It implements an extended version of the reciprocal best alignment heuristic. We apply <monospace>Proteinortho</monospace> to compute orthologous proteins in the complete set of all 717 eubacterial genomes available at NCBI at the beginning of 2009. We identified thirty proteins present in 99% of all bacterial proteomes.</p> <p>Conclusions</p> <p><monospace>Proteinortho</monospace> significantly reduces the required amount of memory for orthology analysis compared to existing tools, allowing such computations to be performed on off-the-shelf hardware.</p

Rapid Pathway Evolution Facilitated by Horizontal Gene Transfers across Prokaryotic Lineages

The evolutionary history of biological pathways is of general interest, especially in this post-genomic era, because it may provide clues for understanding how complex systems encoded on genomes have been organized. To explain how pathways can evolve de novo, some noteworthy models have been proposed. However, direct reconstruction of pathway evolutionary history both on a genomic scale and at the depth of the tree of life has suffered from artificial effects in estimating the gene content of ancestral species. Recently, we developed an algorithm that effectively reconstructs gene-content evolution without these artificial effects, and we applied it to this problem. The carefully reconstructed history, which was based on the metabolic pathways of 160 prokaryotic species, confirmed that pathways have grown beyond the random acquisition of individual genes. Pathway acquisition took place quickly, probably eliminating the difficulty in holding genes during the course of the pathway evolution. This rapid evolution was due to massive horizontal gene transfers as gene groups, some of which were possibly operon transfers, which would convey existing pathways but not be able to generate novel pathways. To this end, we analyzed how these pathways originally appeared and found that the original acquisition of pathways occurred more contemporaneously than expected across different phylogenetic clades. As a possible model to explain this observation, we propose that novel pathway evolution may be facilitated by bidirectional horizontal gene transfers in prokaryotic communities. Such a model would complement existing pathway evolution models

The C-Terminal Domain of the Novel Essential Protein Gcp Is Critical for Interaction with Another Essential Protein YeaZ of Staphylococcus aureus

Previous studies have demonstrated that the novel protein Gcp is essential for the viability of various bacterial species including Staphylococcus aureus; however, the reason why it is required for bacterial growth remains unclear. In order to explore the potential mechanisms of this essentiality, we performed RT-PCR analysis and revealed that the gcp gene (sa1854) was co-transcribed with sa1855, yeaZ (sa1856) and sa1857 genes, indicating these genes are located in the same operon. Furthermore, we demonstrated that Gcp interacts with YeaZ using a yeast two-hybrid (Y2H) system and in vitro pull down assays. To characterize the Gcp-YeaZ interaction, we performed alanine scanning mutagenesis on the residues of C-terminal segment of Gcp. We found that the mutations of the C-terminal Y317-F322 region abolished the interaction of Gcp and YeaZ, and the mutations of the D324-N329 and S332-Y336 regions alleviated Gcp binding to YeaZ. More importantly, we demonstrated that these key regions of Gcp are also necessary for the bacterial survival since these mutated Gcp could not complement the depletion of endogenous Gcp. Taken together, our data suggest that the interaction of Gcp and YeaZ may contribute to the essentiality of Gcp for S. aureus survival. Our findings provide new insights into the potential mechanisms and biological functions of this novel essential protein

Sequence Conservation and Functional Constraint on Intergenic Spacers in Reduced Genomes of the Obligate Symbiont Buchnera

Analyses of genome reduction in obligate bacterial symbionts typically focus on the removal and retention of protein-coding regions, which are subject to ongoing inactivation and deletion. However, these same forces operate on intergenic spacers (IGSs) and affect their contents, maintenance, and rates of evolution. IGSs comprise both non-coding, non-functional regions, including decaying pseudogenes at varying stages of recognizability, as well as functional elements, such as genes for sRNAs and regulatory control elements. The genomes of Buchnera and other small genome symbionts display biased nucleotide compositions and high rates of sequence evolution and contain few recognizable regulatory elements. However, IGS lengths are highly correlated across divergent Buchnera genomes, suggesting the presence of functional elements. To identify functional regions within the IGSs, we sequenced two Buchnera genomes (from aphid species Uroleucon ambrosiae and Acyrthosiphon kondoi) and applied a phylogenetic footprinting approach to alignments of orthologous IGSs from a total of eight Buchnera genomes corresponding to six aphid species. Inclusion of these new genomes allowed comparative analyses at intermediate levels of divergence, enabling the detection of both conserved elements and previously unrecognized pseudogenes. Analyses of these genomes revealed that 232 of 336 IGS alignments over 50 nucleotides in length displayed substantial sequence conservation. Conserved alignment blocks within these IGSs encompassed 88 Shine-Dalgarno sequences, 55 transcriptional terminators, 5 Sigma-32 binding sites, and 12 novel small RNAs. Although pseudogene formation, and thus IGS formation, are ongoing processes in these genomes, a large proportion of intergenic spacers contain functional sequences

Comparative Genomic Analysis of Chitinase and Chitinase-Like Genes in the African Malaria Mosquito (Anopheles gambiae)

Chitinase is an important enzyme responsible for chitin metabolism in a wide range of organisms including bacteria, yeasts and other fungi, nematodes and arthropods. However, current knowledge on chitinolytic enzymes, especially their structures, functions and regulation is very limited. In this study we have identified 20 chitinase and chitinase-like genes in the African malaria mosquito, Anopheles gambiae, through genome-wide searching and transcript profiling. We assigned these genes into eight different chitinase groupings (groups I–VIII). Domain analysis of their predicted proteins showed that all contained at least one catalytic domain. However, only seven (AgCht4, AgCht5-1, AgCht6, AgCht7, AgCht8, AgCht10 and AgCht23) displayed one or more chitin-binding domains. Analyses of stage- and tissue-specific gene expression revealed that most of these genes were expressed in larval stages. However, AgCht8 was mainly expressed in the pupal and adult stages. AgCht2 and AgCht12 were specifically expressed in the foregut, whereas AgCht13 was only expressed in the midgut. The high diversity and complexity of An. gambiae chitinase and chitinase-like genes suggest their diverse functions during different developmental stages and in different tissues of the insect. A comparative genomic analysis of these genes along with those present in Drosophila melanogaster, Tribolium castaneum and several other insect species led to a uniform classification and nomenclature of these genes. Our investigation also provided important information for conducting future studies on the functions of chitinase and chitinase-like genes in this important malaria vector and other species of arthropods

Origin of an Alternative Genetic Code in the Extremely Small and GC–Rich Genome of a Bacterial Symbiont

The genetic code relates nucleotide sequence to amino acid sequence and is shared across all organisms, with the rare exceptions of lineages in which one or a few codons have acquired novel assignments. Recoding of UGA from stop to tryptophan has evolved independently in certain reduced bacterial genomes, including those of the mycoplasmas and some mitochondria. Small genomes typically exhibit low guanine plus cytosine (GC) content, and this bias in base composition has been proposed to drive UGA Stop to Tryptophan (Stop→Trp) recoding. Using a combination of genome sequencing and high-throughput proteomics, we show that an α-Proteobacterial symbiont of cicadas has the unprecedented combination of an extremely small genome (144 kb), a GC–biased base composition (58.4%), and a coding reassignment of UGA Stop→Trp. Although it is not clear why this tiny genome lacks the low GC content typical of other small bacterial genomes, these observations support a role of genome reduction rather than base composition as a driver of codon reassignment